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ABSTRACT 

We  describe  a  preliminary  analysis  of  small-scale  temperature  perturbations  in  the 
stratosphere  produced  by  long  wavelength  gravity  waves,  as  measured  globally  during  the 
CRISTA-SPAS  missions.  We  focus  on  enhanced  activity  observed  at  the  equator  and  near  the 
southern  tip  of  South  America.  Global  ray-tracing  simulations  for  the  mission  days  indicate  that 
the  equatorial  activity  is  broadly  consistent  with  the  transmission  of  nonzero  phase  speed  waves 
into  the  stratosphere.  Ray-tracing  and  mountain  wave  model  simulations  indicate  that  the 
activity  over  South  America  is  produced  by  long  wavelength  mountain  waves  forced  by  flow 
over  the  Andes.  The  results  suggest  that  study  of  global  gravity  wave  “weather”  in  the  middle 
atmosphere  may  soon  be  feasible. 

1.  INTRODUCTION 

A  range  of  weather-related  phenomena  in  the  lower  atmosphere  generates  gravity  waves. 
Examples  include  fronts  (e.g.,  Eckennann  and  Vincent  1993),  squall  lines  (Alexander  et  al. 
1995),  convective  clouds  (e.g.,  Alexander  and  Pfister  1995),  cyclogenesis  (e.g.,  Powers  1997), 
typhoons  (e.g.,  Sato  1993),  geostrophic  adjustment  of  jet  steams  (e.g.,  Kaplan  et  al.  1997),  flow 
over  mountains  (e.g.,  Ralph  et  al.  1993),  and  so  on.  If  wind  patterns  aloft  are  favorable,  some  of 
these  waves  can  propagate  into  the  middle  atmosphere.  Due  to  decreasing  atmospheric  densities, 
gravity  waves  grow  in  amplitude  with  height  and  so  must  eventually  generate  instabilities  and 
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“break.”  Such  quasi-continuous  gravity  wave  breaking  generates  drag  and  diffusion  that 
strongly  influence  the  global  circulation  and  temperature  patterns  of  the  middle  atmosphere 
(e.g.,  Haynes  et  al.  1991;  Fritts  and  Luo  1995;  Alexander  and  Holton  1997). 

While  the  impact  of  breaking  gravity  waves  on  the  climatology  of  the  middle 
atmosphere  is  well  appreciated,  far  less  is  known  about  how  shorter-term  changes  in  weather- 
related  gravity  wave  sources  and  local  propagation  environments  affect  the  “weather”  in  the 
middle  atmosphere.  Such  research  has  been  stymied  to  date  by  a  lack  of  data  on  the  synoptic 
meteorology  of  the  middle  atmosphere  and  its  contemporaneous  gravity-wave  activity.  There 
have  also  been  very  few  global  models  of  gravity-wave  production  and  propagation  with  any 
sort  of  forecasting  capability. 

These  shortcomings  in  observations  and  modeling  are  slowly  being  resolved,  allowing 
us  to  investigate  “gravity  wave  weather”  issues  systematically  for  the  first  time.  Next  generation 
satellite  instruments  have  sufficient  spatial  resolution  to  resolve  explicitly  large-scale  gravity 
waves  in  the  middle  atmosphere  (e.g.,  Ross  et  al.  1992;  Fetzer  and  Gille  1994;  Mende  et  al. 
1994;  Wu  and  Waters  1996;  Dewan  et  al.  1998;  Preusse  et  al.  1999;  Eckermann  and  Preusse 
1999).  Similar  advances  in  modeling  have  led  to  global  models  that  can  simulate  and  even 
forecast  gravity  wave  activity  generated  by  specific  sources  (e.g.,  Bacmeister  et  al.  1994;  Marks 
and  Eckermann  1995;  Alexander  and  Holton  1997;  Dombrack  et  al.  1998;  Sato  et  al.  1999; 
Charron  and  Brunet  1999).  This  has  enabled  some  initial  comparisons  between  global  models 
and  global  measurements  of  gravity  waves  in  the  middle  atmosphere  (Eckennann  and  Marks 
1997;  Alexander  1998;  Eckennann  and  Preusse  1999). 

We  focus  here  on  preliminary  data  from  two  observational  campaigns  with  the  CRISTA- 
SPAS  satellite  during  shuttle  missions  STS-66  and  STS-85  (November,  1994,  and  August, 
1997,  respectively).  Analysis  shows  that  large-scale  gravity  wave  fluctuations  are  present  in  the 
stratospheric  and  mesospheric  temperature  measurements.  These  conclusions  are  buttressed  by 
numerical  models,  which  generate  maps  of  gravity  wave  activity  that  resemble  the  measured 
fluctuations.  These  results  provide  some  preliminary  insights  into  the  generation  and  filtering 
processes  affecting  gravity  waves  in  the  middle  atmosphere. 

2.  THE  CRISTA-SPAS  MISSIONS 

The  Cryogenic  Infrared  Spectrometers  and  Telescopes  for  the  Atmosphere  (CRISTA) 
instrument  operated  aboard  the  Shuttle  Pallet  Satellite  (SPAS)  (Offermann  et  al.  1999). 
CRISTA-SPAS  was  deployed  into  orbit  and  later  recovered  by  the  shuttle  during  missions  STS- 
66  (November  3-14,  1994)  and  STS-85  (August  7-16,  1997).  Figure  1  shows  deployment  of 
CRISTA-SPAS  by  the  orbiter’s  robot  ann  on  November  4,  1994.  Lee  wave  clouds  downstream 
(west)  of  the  Andean  ridge  over  the  Chubut  River  region  in  central  Patagonia  (Argentina)  are 
visible  in  this  photograph,  indicating  that  mountain  waves  were  being  generated  by  flow  over 
the  Andes  as  CRISTA-SPAS  was  being  prepared  for  orbit. 

CRISTA  is  a  limb-scanning  instrument  that  measures  infrared  emission  spectra  with 
high  spectral  resolution  (Riese  et  al.  1999;  Grossmann  2000).  Temperatures  in  the  stratosphere 
and  lower  mesosphere  are  retrieved  from  the  CCL  Q-branch  emission  at  792  cm'1  (-12.6  pm). 
Full  details  of  the  “version  1”  retrievals  are  given  by  Riese  et  al.  (1999),  which  yielded 
temperatures  with  a  precision  -1  K.  The  precision  of  the  “version  2”  data  analyzed  here  is  ~0.5 
K.  CRISTA  uses  two  additional  telescopes  that  look  18°  either  side  of  the  central  telescope, 
corresponding  to  tangent  heights  separated  by  -650  km.  Each  scanning  telescope  acquires  a 
complete  vertical  profile  every  200-400  km  along  track,  depending  on  the  measuring  mode. 
These  unique  measurements  yielded  data  with  global  coverage  and  high  spatial  resolution  (see 


Figure  2  of  Riese  et  al.  1999). 

Mesospheric  temperatures  are 
retrieved  from  CO2  emissions 
in  the  652-657  cm'1  band, 
which  are  measured  by  the 
center  telescope  only.  Their 
precision  is  estimated  to  be 
~1  K.  The  vertical  sampling 
interval  is  1.5  km  for  both  the 
stratospheric  and 

mesospheric  temperature 
data.  Further  details  of  the 
instrument  and  measurements 
are  given  by  Offermann  et  al. 

(1999),  Riese  et  al.  (1999) 
and  Grossmann  (2000). 

Fetzer  and  Gille 
(1994,  1996)  found  large- 
scale  gravity-wave 

fluctuations  in  stratospheric 
data  acquired  by  the  Limb 
Infrared  Monitor  of  the 
Stratosphere  (LIMS),  which 
operated  on  the  Nimbus  7 
satellite  during  1978-79.  The 
precision  and  resolution  of 
CRISTA  measurements 
exceed  those  of  LIMS. 

Furthermore,  CRISTA  measures  spectra  whereas  LIMS  was  a  filter  instrument.  This  allows 
CRISTA  to  check  the  quality  of  the  retrievals  in  regions  of  enhanced  aerosol  content  (e.g.,  the 
tropical  lower  stratosphere).  Retrieval  modeling  by  Preusse  et  al.  (1999)  indicates  that 
sinusoidal  gravity  wave  temperature  perturbations  with  wavelengths  in  excess  of  -5  km 
vertically  and  -200  km  horizontally  are  resolved  by  the  CRISTA  temperature  channels.  Thus, 
gravity  waves  with  wavelengths  longer  than  these  limits  should  be  evident  in  the  CRISTA 
temperature  data. 

To  investigate  this,  we  isolated  small-scale  temperature  variability  in  the  CRISTA 
temperature  data.  Synoptic  variability  was  removed  using  a  wavenumber  0-6  Kalman  filter. 
Next,  vertical  profiles  of  residual  fluctuations  were  analyzed  using  the  Maximum  Entropy 
Method  (MEM)  and  harmonic  analysis.  MEM  analysis  was  used  to  identify  the  energetic 
vertical  wavelength  oscillations  in  a  given  profile.  These  wavelengths  were  then  used  to 
constrain  a  harmonic  fit  over  a  12  km  altitude  interval,  which  was  progressively  shifted  over  the 
full  height  range  of  the  profile.  In  this  way,  two  principal  amplitudes  and  phases  were 
determined  as  a  function  of  height  for  each  profile.  Finally,  a  multiplier,  determined  by  the 
instrument’s  “visibility”  to  this  particular  vertical  wavelength,  was  used  to  scale  the  fitted 
amplitude.  For  further  details,  see  Preusse  et  al.  (1999)  and  Eckermann  and  Preusse  (1999). 

Figure  2  shows  resulting  global  maps  of  small-scale  temperature  amplitudes  at  20-30  km 
for  6  successive  days  during  Mission  1  (STS-66).  We  note  a  band  of  enhanced  variance  at  the 
equator  that  exhibits  significant  zonal  asymmetry  and  day-to-day  variability.  This  is 


Figure  1:  Astronaut  photo  of  CRISTA-SPAS  being  deployed  into  orbit  by 
the  robot  arm  of  the  Space  Shuttle  Atlantis  at  12:19  UT  on  4th  November, 
1994.  Visible  below  are  lee  wave  clouds  downstream  of  the  southern 
Andes  over  the  Chubut  region  of  Argentina  (~42°S).  This  photograph  and 
others  from  STS-66  are  accessible  from  NASA’s  shuttle  image  archive,  at 
http://images.jsc.nasa.gov/iams/images/earth/STS066/html/20177512.htm. 
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Figure  2:  Preliminary  global  maps  of  CRISTA  temperature  amplitudes  in  the  20-30  km  height  range  for  six 
successive  days  during  Mission  1  (5th- 10th  Nov 
and  blue-white  shading  shows  smallest  values. 


successive  days  during  Mission  1  (5th- 10th  November,  1994).  Red  (dark)  shading  indicates  highest  values  (~2.5  K), 


superficially  consistent  with  the  notion  of  enhanced  convective  launching  of  gravity  waves  and 
equatorial  waves  into  the  tropical  stratosphere  (e.g.,  Alexander  and  Pfister  1995;  Garcia  2000). 
We  also  note  a  persistent  zone  of  enhanced  activity  near  the  southern  tip  of  South  America, 
consistent  with  the  concept  of  mountain  waves  generated  by  flow  over  the  southern  Andes  (see 
Figure  1)  penetrating  into  the  stratosphere.  However,  work  by  Alexander  (1998)  has 
demonstrated  that  straightforward  interpretations  of  raw  satellite  gravity  wave  data  such  as  these 
can  be  potentially  misleading.  Therefore  we  turn  to  global  numerical  modeling  of  gravity  waves 
to  make  further  progress  in  interpreting  these  and  other  features. 


3.  GLOBAL  GRAVITY  WAVE  MODELS 


3.1.  Gravity  Wave  Regional  or  Global  Ray  Tracer  (GROGRAT) 

GROGRAT  is  a  numerical  model  that  uses  ray-tracing  techniques  to  track  the 
propagation  and  amplitude  evolution  of  nonhydrostatic  gravity  waves  within  arbitrary  regional 
or  global  specifications  of  the  Earth’s  atmosphere.  Marks  and  Eckermann  (1995)  described  the 
basic  formulation  of  the  model,  while  Eckennann  and  Marks  (1997)  reported  subsequent 


improvements  and  updates  to  the  code.  Further  information  can  be  found  at  the  GROGRAT  web 
site1. 

For  the  background  atmosphere,  we  use  6  hourly  global  2.5°x2°  (144x91)  wind  and 
temperature  data  from  NASA’s  Data  Assimilation  Office  (DAO),  which  come  gridded  at  18 
standard  pressure  levels  from  1000  to  0.4  hPa  (Schubert  et  al.  1993;  Coy  and  Swinbank  1997). 
We  resampled  these  data  onto  a  64x45x18  global  grid.  Comparison  of  Figures  3a  and  3b  shows 
that  the  major  weather  phenomena  are  retained  in  the  wind  fields  after  this  resampling. 
GROGRAT  absorbed  these  fields,  first  regridding  them  onto  a  geometrical  height  grid  (using 
the  DAO  geopotential  height  fields),  then  fitting  the  data  at  each  height  with  a  superposition  of 
spherical  harmonics  Ym  n((p,0).  Lower  panels  in  Figure  3  show  fits  using  spherical  harmonics  out 
to  orders  (m,n)=(  16,16)  and  (m,n)=( 32,32)  (panels  c  and  d,  respectively).  Both  capture  the 
overall  structure  quite  well,  but  Figure  3d  captures  more  of  the  finer-scale  weather-related 
details  in  the  winds.  Thus,  we  used  fits  out  to  {m,n)={ 32,32)  in  our  simulations.  Note  also  that 
GROGRAT  interpolates  all  the  individual  spherical  hannonic  coefficients  at  each  fitting  level 
with  cubic  splines,  which  ensures  smoothly  varying  continuous  derivatives  in  all  three  spatial 
directions. 

Here  we  ignore  sources  and 
concentrate  on  the  transmission 
characteristics  of  the  atmosphere.  To 
this  end,  we  used  GROGRAT  to  trace 
a  globally  unifonn  distribution  of 
waves  through  the  middle  atmosphere 
using  DAO  winds  and  temperatures 
during  the  mission  days.  Waves  were 
launched  at  zo  =1  km  with  ground- 
based  phase  speeds  c  =  0,  20  and  40  m 
s'1,  a  horizontal  wavelength  Ax  =  50 
km,  a  peak  initial  horizontal  velocity 
amplitude  of  u  ’  =  0.2  m  s'1  and  an 
isotropic  distribution  of  8  initial 
propagation  directions.  Wave 
amplitudes  were  controlled  using  a 
wave  action  conservation  equation 
with  damping  terms  due  to  turbulent 
diffusion,  infrared  radiative  cooling 
and  wave  breaking/saturation  (Marks 
and  Eckermann  1995).  A  total  of 
13056  rays  were  traced  over  the  globe. 

Figure  4  shows  global  ray 
counts,  weighted  by  the  horizontal 
velocity  amplitude  of  each  ray,  as  a 
function  of  height  on  6th  November, 
1994  at  00:00  UT.  Results  shown  are 
for  the  c  =  20  m  s'1  waves,  as  might 
plausibly  emanate  from  convection  and 
other  moving  and/or  evolving  weather 


GSFC  Data  Assimilation  Office  (DAO)  Analyses:  Nov  5 1994 
Zonal  Winds  at  48°S 


(longitude-height) 

(a)  DAO  Original:  144x91x18  (b)  DAO  Resampled:  64x45x18 


Figure  3:  Longitude-height  plot  of  DAO  zonal  wind  at  48°S  on 
5th  November,  1994  at  00:00  UT:  (a)  original  data;  (b) 
resampled  to  64x45  in  the  horizontal;  (c)  after  GROGRAT 
regridding  to  geometric  heights  and  spherical  harmonic  fits  to 
maximum  (m,n)=(  1 6, 1 6);  (d)  as  for  (c),  but  after  spherical 
harmonic  fits  to  maximum  (w,7?)=(32,32).  Contour  labels  are  in 
m  s"1,  positive  values  are  eastward  and  have  solid  contour  lines. 


1  http://uap-www.nrl.navy.mil/dynamics/html/grograt.html 


Figure  4:  GROGRAT  simulations  of  global  gravity  wave  activity  on  6th  November,  1994,  using  a  globally 
invariant  source  spectrum  of  waves  of  c=20  m  s"1  (see  text).  Light  values  indicate  large  wave  variances,  darker 
values  indicate  small  wave  variances. 


systems.  Although  wave  amplitudes  are  globally  uniform  at  z0  =1  km,  the  variable  background 
atmosphere  quickly  yields  considerable  geographical  variability  at  upper  heights.  At  z  =  20  km 
and  above,  Figure  4  reveals  a  clear  band  of  enhanced  activity  around  the  equator  that  resembles 
the  enhanced  equatorial  activity  at  20-30  km  in  the  CRISTA  data  (Figure  2).  It  occurs  in  the 
model  due  to  very  light  winds  throughout  the  equatorial  troposphere  and  stratosphere  during 
Mission  1,  which  enables  waves  with  nonzero  phase  speeds  to  propagate  upwards  without 
encountering  critical  layers  or  severe  damping.  However,  the  zonal  asymmetry  in  Figure  2  is  not 
reproduced  in  Figure  4,  suggesting  that  zonally  varying  wave  generation  may  also  play  a  role  in 
producing  the  observed  activity. 

We  are  investigating  possible  sources  of  equatorial  wave  activity  during  the  missions. 
For  example,  during  Mission  2  (STS-85;  7-16  August,  1997)  similar  bands  of  enhanced 
equatorial  activity  were  observed.  Super  Typhoon  Winnie  developed  over  the  South  China  Sea 
during  the  mission,  reaching  Category  5  (“super  typhoon”  status)  on  August  12  (see  Figure  5). 
Typhoons  are  known  to  be  strong  sources  of  gravity  waves  (e.g.,  Matsumoto  and  Okamura 
1985;  Sato  1993),  and  thus  Winnie  might  have  produced  enhanced  stratospheric  gravity  wave 
activity.  CRISTA  used  a  “hawkeye”  observational  mode  during  Mission  2,  which  increased 
sampling  over  Indonesia  but  reduced  overall  coverage  over  the  typhoon.  Nevertheless, 
preliminary  observations  in  this  region  show  some  suggestions  of  enhanced  temperature 


Figure  5:  (a)  full  disk  GOES  8  visible  image  on  12m  August,  1997,  showing  Super  Typhoon  Winnie.  Australia  is 
visible  to  the  south  (image  courtesy  ofNOAA’s  National  Climatic  Data  Center);  (b)  close-up  photo  of  Winnie 
taken  from  Discovery  on  13th  August  (robot  arm  in  foreground). 


fluctuations  in  the  stratosphere  near  the  typhoon.  Investigation  of  equatorial  source  processes 
continues. 


3.2.  Naval  Research  Laboratory  Mountain  Wave  Forecast  Model  (MWFM) 

Mountains  are  believed  to  be  a  strong  extratropical  source  of  gravity  waves  for  the 
middle  atmosphere  (Nastrom  and  Fritts  1992;  Bacmeister  1993).  With  this  in  mind,  the  Naval 
Research  Laboratory  Mountain  Wave  Forecast  Model  (MWFM)  was  developed  to  forecast  the 
turbulence  produced  by  breaking  mountain  waves  in  the  stratosphere,  as  an  aid  to  safe  flight 
planning  for  NASA’s  stratospheric  ER-2  research  aircraft  (Bacmeister  et  al.  1994;  Eckermann  et 
al.  2000).  The  model  uses  a  database  of  quasi-two-dimensional  ridges  inferred  from  high- 
resolution  digital  topography.  Surface  winds  from  forecast  data  are  blown  across  these  ridges 
and  used  to  generate  waves.  Wave  equations  are  used  thereafter  to  trace  any  vertical  extension 
of  this  activity  into  the  atmosphere,  as  well  as  any  subsequent  wave  breaking.  Further 
information  can  be  found  at  the  MWFM  web  site  . 

Here  we  use  the  MWFM  with  DAO  data  to  generate  an  initial  global  field  of  mountain 
waves  forced  at  the  ground.  We  use  GROGRAT  thereafter  to  track  the  propagation  of  these 
waves  to  greater  heights,  using  fitted  DAO  data  as  before.  Figure  6  shows  modeled  mountain 
wave  activity  at  z  =  25  km  from  5-10  November,  1994.  We  note  that  stratospheric  mountain 
wave  activity  is  predicted  at  the  southern  tip  of  South  America  during  Mission  1,  which 
resembles  the  activity  observed  by  CRISTA  in  Figure  2.  However,  there  is  additional  activity 
over  Eurasia  and  western  North  America  in  Figure  6  that  is  not  seen  as  clearly  in  Figure  2. 

In  analyzing  MLS  radiance  fluctuations,  Alexander  (1998)  stressed  the  importance  of 
accounting  for  the  spatial  filtering  of  the  satellite  gravity  wave  measurement.  Mountain  wave 
vertical  wavelengths  are  given  approximately  by 


2  http://uap-www.nrl.navy.mil/dynamics/html/mwfm.html 


(1) 


X,  = 


2nU  cos  (p 
N 


where  U  is  the  horizontal  wind  speed,  N  is  the  Brunt- Vaisala  frequency  and  cp  is  the  angle 
between  the  horizontal  wind  and  wave  vectors.  Between  40°-50°S  over  South  America  (— 70°W) 
the  DAO  winds  blow  approximately  orthogonal  to  the  Andean  ridge  axis,  implying  cp  -  0°.  At 
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Figure  6:  Coupled  GROGRAT-MWFM  simulations  of  global  mountain  wave  amplitudes  at  z  =  25  km,  using 
DAO  winds  and  temperatures.  Light  (dark)  values  depict  large  (small)  wave  activity. 


20-30  km,  U  -20-50  m  s’1  (see  Figure  3)  and  N  ~  0.02  rad  s’1,  yielding  -  6-16  km  from  (1). 
Conversely,  U  -10  m  s’1  at  30°N,  yielding  vertical  wavelengths  no  greater  than  -3  km  here. 
Since  CRISTA  is  sensitive  to  vertical  wavelengths  Xz>  5  km,  it  can  resolve  the  long  Xz  waves 
over  South  America  but  not  the  shorter  Xz  waves  at  northern  latitudes. 

To  assess  this  quantitatively,  Figure  7  shows  the  model  results  from  Figure  6  after 
imposition  of  a  filter  which  retains  only  those  waves  with  vertical  wavelengths  Xz  >  10  km.  We 
see  that  the  activity  over  southern  South  America  is  now  more  prominent  and  activity  in  the 
Northern  Hemisphere  is  significantly  attenuated,  which  is  more  in  line  with  the  data  in  Figure  2. 
Alexander  (1998)  reported  similarly  improvements  on  comparing  filtered  model  data  to  the 
MLS  radiance  fluctuations  reported  by  Wu  and  Waters  (1996). 


Amplitude-Weighted  Ray  Counts:  z=25  km 


Figure  7:  As  in  Figure  6,  but  after  imposition  of  a  filter  which  retains  A>  10  km  waves  only. 


Eckermann  and  Preusse  (1999)  present  an  in-depth  study  of  temperature  fluctuations 
observed  by  CRISTA  over  southern  South  America  and  central  Eurasia.  Using  simple  mountain 
wave  theory  and  “next  generation”  MWFM  hindcast  simulations,  they  show  conclusively  that 
this  activity  is  produced  by  long  wavelength  mountain  waves  which  propagate  into  the 
stratosphere.  Mesoscale  modeling  work  by  Tan  and  Eckermann  (2000)  generates  stratospheric 
mountain  waves  very  similar  to  those  measured  by  CRISTA  over  the  southern  Andes  on  6th. 
November,  1994. 

4.  SUMMARY 

Initial  comparisons  between  the  small-scale  temperature  perturbations  measured  by 
CRISTA  (Figure  2)  and  global  gravity  wave  model  predictions  are  encouragingly  favorable. 
These  findings  further  support  the  conclusion  that  CRISTA  temperature  measurements  in  the 
stratosphere  and  mesosphere  resolved  large-scale  gravity  waves  (see  also  Preusse  et  al.  1999; 
Eckermann  and  Preusse  1999).  Global  synoptic  data  on  gravity  waves  in  the  middle  atmosphere 
are  still  very  limited,  having  been  inferred  from  a  few  earlier  satellite  measurements  only 
(Fetzer  and  Gille  1994  1996;  Wu  and  Waters  1996).  Furthermore,  interpretation  of  some  of 
these  earlier  measurements  has  not  always  been  straightforward  (Alexander  1998).  Yet  global 
measurements  of  gravity  waves  and  accompanying  theoretical  explanations  of  the  observed 
features  are  vital  to  a  proper  understanding  of  the  “weather”  in  the  middle  atmosphere,  which  is 
very  sensitive  to  gravity  wave  driving  from  below.  Our  analysis  suggests  that  some  key  features 


in  the  CRISTA  data  are  consistent  with  specific  gravity  wave  production  and  filtering  processes. 

Further  analysis  and  modeling  of  these  data  are  being  pursued. 
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